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Abstract 



In a class of recently proposed models, the early universe is strongly coupled and described 
holographically by a three-dimensional, weakly coupled, super-renormalizable quantum field 
theory. This scenario leads to a power spectrum of scalar perturbations that differs from the 
usual empirical ACDM form and the predictions of generic models of single field, slow roll 
infiation. This spectrum is characterized by two parameters: an amplitude, and a param- 
eter g related to the coupling constant of the dual theory. We estimate these parameters, 
using WMAP and other astrophysical data. We compute Bayesian evidence for both the 
holographic model and standard ACDM and find that their difference is not significant, al- 
though ACDM provides a somewhat better fit to the data. However, it appears that Planck 
will permit a definitive test of this holographic scenario. 
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1. Introduction 

It was recently proposed that the very early universe is in a strongly coupled and non- 
geometric phase, best described holographically as a weakly coupled three-dimensional quan- 
tum field theory in flat spacetime with no gravity [T||2||3|H] . In this scenario, the strongly 
coupled phase plays the same role as inflation in more conventional models of the very early 
universe. The corresponding predictions for the scalar and tensor power spectra and non- 
Gaussianity were worked out in [2] and [1] (see [3] for a summary). In the present paper, 
we confront these models with the observational data from WMAP and other sources, and 
compare them to the fit provided by the standard power-law ACDM model. 

This holographic scenario is qualitatively different from conventional models of inflation. 
The latter are described using graviton and inflaton fluctuations about a homogeneous and 
isotropic inflating background. In the new holographic models, however, the very early 
universe is such that notions of spacetime, and perturbations around it, are not yet well- 
defined. Rather, they emerge at the end of the strongly coupled epoch, after which the 
universe is described by usual hot big bang cosmology. Such a non-geometric period should 
have a well-defined description in string theory in terms of a strongly coupled sigma model; 
here, we use holography to model it in terms of a weakly coupled three-dimensional QFT. 
Conventional inflationary models can also be described using the holographic framework, 
but these scenarios are dual to strongly coupled QFTs. 

A wealth of observational data now allows us to test theoretical ideas about the early 
universe, providing a probe of fundamental physics at energies close to the Planck scale. 
Power-law ACDM, an empirical model with only six free parameters, fits the existing data 
remarkably wellll] This model describes a flat universe with radiation, baryons, cold dark 

^Cosmological fits to CMB data usually include a parameter quantifying the contribution of the Sunyaev- 
Zel'dovich effect to the apparent small-scale temperature anisotropies. We marginalize over this parameter, 
but will not discuss it further in what follows. 
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matter, a cosmological constant and a power-law spectrum of adiabatic primordial fluctua- 
tions. Four of the six parameters describe the composition and expansion of the universe, 
namely the Hubble rate Hq = lOO/ikm/s/Mpc, the physical baryon and dark matter densi- 
ties fife/i^ and flch'^, and the optical depth due to re- ionization r. Given that we do not need 
spatial curvature to fit the data, the current dark energy contribution follows from the re- 
quirement that the overall density of the universe is equal to the critical value. The remaining 
two parameters characterize the power spectrum of primordial curvature perturbations 



where A^(g^,) is the amplitude and Ug the spectral tilt. The measured value of A^(g*) 
depends on an arbitrary reference scale, or pivot, q^,. 

This is a purely empirical parameterization but there is no evidence that the primordial 
power spectrum is not well-described by this choice (see e.g. |5j). Moreover, while almost all 
slow roll inflationary models predict that Ug is a function of q, the dependence is typically 
weak. This "running" is expressed in terms of as = dus/dlnq, evaluated at g*. In simple 
inflationary models, is higher order in slow roll than — 1, the departure from scale 
invariance [6]. 

Interestingly, the holographic model predicts a scalar power spectrum of the form 



where g is & free parameter that replaces Ug. This power spectrum corresponds to a 2- loop 
approximation on the QFT side, and higher loop terms would modify this formula. Roughly 
speaking, gq^^/q corresponds to the effective coupling constant, and for self-consistency it 
should be small for all cosmologically relevant values of q. When this is the case, one can 
transform (II. 2p into the form of (11. ip with a spectral index rig = ns{q) that is now scale- 
dependent [2] . In this case, however, all logarithmic derivatives d^Ug (q) /d\nq^ are of the same 
order, in contrast to the slow roll case, and the expansion cannot be truncated. A principal 
goal of this paper is to determine whether a spectrum of the form (II. 2p is consistent with 
the data and to quantify its ability to explain the observed universe, relative to that of the 
conventional power spectrum. In doing so, we are contrasting the predictions of holographic 
cosmology with both the empirical ACDM model and the large class of slow roll inflationary 
scenarios whose power spectrum is effectively described by (II. ID . The parameter set for 
"holographic ACDM" is (ytbh"^, ^ch"^, h, r, A^, g), while the ACDM parameter set takes its 
usual form. 

The form of (II. 2p follows uniquely from the basic properties of a certain class of dual 
QFTs. If the data are incompatible with (11.21) we could then exclude this entire class of holo- 
graphic models. More precisely, the dual model is a super-renormalizable QFT admitting a 
't Hooft large- iV limit, with massless fields and a single dimensionful coupling constant. The 
momentum dependence of (11.21) essentially follows from dimensional analysis: the effective 
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coupling of the theory is obtained by forming the unique dimensionless combination of the di- 
mensionful couphng constant and the momentum q. As mentioned above, (11. 2p corresponds 
to a 2-loop approximation, and all momentum dependence enters solely through the depen- 
dence of this computation on the effective coupling. In particular, the scale-invariance of 
(II. 2p at large q is directly related to the asymptotic freedom of super-renormalizable theories. 

Given that the computation of (II. 2p rests on two approximations, the large- limit and 
a perturbative 2-loop approximation, in addition to asking whether holographic ACDM fits 
the data, we must also check that the best-fit values of and g are compatible a posteriori 
with these theoretical approximations. The smallness of the amplitude immediately 
implies that A^ ~ 10^, so the large-A^ approximation is indeed valid. Moreover, the best-fit 
value of g is very small, reflecting the near scale-invariance of the power spectrum, but it will 
be important to check that the effective coupling is small for all values of q which contribute 
to the observed CMB anisotropics and are thus constrained by the WMAP data. 

Our data analysis is performed using a suitably modified version of CAMB and Cos- 
moMC [7ll8] . together with the MultiNest sampler |9llT0] . As well as speeding the parameter 
estimations relative to the usual MCMC sampler^, this combination of tools efficiently com- 
putes Bayesian evidence. Evidence discriminates between two or more competing models, 
and is sensitive to the likelihood across the overall parameter space of each model p!T|T2] . 
In addition, we compare the maximum likelihood values found for each model, and the 
corresponding information criteria. 

This paper is organized as follows. In Section 2, we summarize the theoretical calculation 
of the power spectrum for holographic ACDM. Section 3 examines the parameter estimates 
obtained for this model and the appropriate priors for the free parameters. Section 4 discusses 
the model selection problem, and Section 5 contains a discussion of our results. 

As this paper was being completed [13] appeared, which also presents parameter con- 
straints for the holographic ACDM power spectrum. However, the parameter estimation 
problem is approached somewhat differently here, and we also address model selection. 

2. Theoretical prediction 

In this section we describe the theoretical prediction for the power spectrum of the holo- 
graphic model. 

In general, a holographic model is specified by providing the dual three-dimensional 
QFT and the holographic formulae that relate the cosmological observables of interest to 
correlation functions of the dual QFT. There are two classes of three-dimensional theories 
that one may currently use to model the very early universe: the first class consists of three- 
dimensional QFTs with a non-trivial fixed point in the UV, while the second class involves 

^The performance of an MCMC sampler is a function of the allowed parameter ranges, stepsizes, data 
and the underlying model. We could not find a set of parameters for the usual MCMC sampler with the 
holographic model for which the acceptance rate was close to that typically found for ACDM, resulting in 
chains that converged only very slowly. 
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Figure 1: 1-loop contribution to 
{T{q)T{—q)). We sum over the 
contributions from gauge fields, 
scalars and fermions, witfi eacfi 
diagram yielding a contribution 
of order ~ N'^q^. 



super-renormalizable QFTs. The models we discuss here belong to this second class. We 
emphasize that the theoretical predictions discussed below are specific to this class of models. 

More precisely then, a theory in this second class should have the following properties: 
(i) it should admit a large- limit, (ii) all fields should be massless, (iii) it should have 
a dimensionful coupling constant, (iv) all terms in the Lagrangian should have the same 
scaling dimension, which should be different from three. Properties (ii)-(iv) imply that the 
theory admits a generalized conformal structure p^pS] , meaning that the theory would 
be conformal if the coupling constant is promoted to a background field that transforms 
under conformal transformations. A class of models exhibiting these properties is given by 
three-dimensional SU{N) Yang-Mills theory coupled to a number of massless scalars and 
fermions, all transforming in the adjoint of SU{N), with interactions consisting of Yukawa 
terms and quartic scalar terms. In three dimensions, the Yang-Mills coupling constant Qym 
has dimension one, and for the models we discuss, one may arrange (by rescaling the fields 
appropriately) that the coupling constant appears only as an overall constant in the action. 
Assigning scaling dimension one to scalars and gauge fields, and 3/2 to fermions, one finds 
that kinetic terms and the interactions have dimension four. Theories of this type appear as 
the worldvolume theories of D-branes. 

In our case, we are interested in the primordial power spectrum of scalar perturbations: 
this may be obtained from the 2-point function for the trace of the stress tensor T as we 
discuss below. The generalized conformal structure and large-A^ counting imply that the 
general form of the 2-point function at large A^ i^ [15], 



(T(g)T(-g)) = q'N'f{g!^), (2.1) 

where q is the magnitude of the 3-momentum, g^^ = g\y[N/q is the effective dimensionless 
't Hooft coupling and fig^s) is a general function of g^^. 

The holographic framework of [1] involves the analytic continuation 

q -iq, N -iN, (2.2) 

and the formula for scalar power spectrum is 

^niQ) = -i^im(T(g)r(-g))' ^^'^^ 



■^If one imposes only the generalized conformal structure then the r.h.s. of (|2.ip is modified as N'^fig'^g) — > 
f{N^, ffoff)' where f{N^, g^g) is a general function of two variables. 
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Figure 2: Diagram topologies contributing at 2-loop order. 



where the imaginary part is taken after applying the analytic continuation. Note that under 
(ESD we have 

^7eff -> als. N'q' ^ -tN'q\ (2.4) 

It follows that in theories with generalized conformal invariance, the 2-point function of the 
stress tensor transforms very simply under the analytic continuation (12. 2p . We thus obtain 
the final formula for the scalar power spectrum: 

^Uq) = T-^TTir- (2.5) 



2 

cff^ 



Including subleading l/N"^ corrections to this formula is straightforward (one just expands 
the function /(A^^, g^f^) in l/N^, see footnote|3]), but, as we shall see, the data favor ~ 10^ 
so such terms are in practice neglig ibleE 

When (7gfj is small, one may compute the function f{g^g) and one finds that it has the 
form 

/(^?eff) = /0(1 - /l^7eff Inge's + /2^7eff + 0[gts]). (2.6) 

The leading-order term /q is determined at 1-loop in perturbation theory (see Fig. [1]), and 
its precise value may be found in |2]. The constant /i is determined by a standard 2-loop 
computation. As is well known, in perturbation theory super-renormalizable theories with 
massless fields display severe infrared divergences. Indeed, each of the 2-loop diagrams listed 
in Fig. [2] evaluates to an overall factor of A^^^fyM multiplying an integral with superficial 
degree of (infrared) divergence two. Imposing an infrared cut-off, giR, one may evaluate the 
integrals to obtain ~ g^ln(g/giR). Altogether, one finds a 2-loop contribution to the stress 
tensor 2-point function of the order 



N'q'g'^s Hq/Qir) = N'q\-gi^ \ngi^ + gi^ \n{g^^N / qi^)) . (2.7) 

Thus, /i is determined but /2 is still undetermined since so far giR is arbitrary. It was argued 
in [16], however, that this infrared divergence is an artefact of perturbation theory and the 
theory develops a physical scale that acts as a cut-off. In [T7], it was shown for a specific 



^ On the other hand, it is unclear how to incorporate effects that are non-perturbative in N: exponen- 
tiahy smah corrections of the type exp(— A^^) are mapped to exponentiahy large effects after the analytic 
continuation. In the discussion above it is understood that the large- limit is taken before the analytic 
continuation is performed. 
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Figure 3: Perturbative theoretical prediction for the power spectrum of the holographic 
model. The lower curve corresponds to g > Q while the upper corresponds to g < 0. The 
perturbative calculation is rehable for gl^ ~ g<i*/<i ^ 1, corresponding to large momenta 
i/qI* ^ 1 far from the peak/trough feature at In \q/gq*\ = 1. At sufficiently high momenta, 
the power spectrum becomes nearly scale invariant, with g > corresponding to a blue tilt 
and < to a red tilt. 

class of models that a large- resummation indeed leads to a finite answer with giR ~ g\y^N . 
Having thus obtained the infrared scale, /2 is then determined unambiguously. 
By rearrangement, (12 .6^ may also be written as 



where /s = exp(— /2//1). The constants /i, /2, /s in general depend on all parameters of the 
theory {i.e., the field content, Yukawa and quartic couplings), but have not been computed 
to date. As long as we probe the theory at scales far above the infrared scale giR, however, 
the specific value of /s should only provide a small correction since | In g^^ \ ^ | In /s | . We 
will thus write /s = and take /3 = 1 in the following; later, in Section H73l we will check 
how the fit to the data changes if (3 is allowed to vary. To simplify our notation, we set 



where g* is the pivot scale, taken here to be g^, = 0.05 Mpc ^. Substituting back into 
(12. Sp . we finally obtain the power spectrum written down in equation (II. 2p after making 
the identification = 1 / {Att^ fo) . Note that in previous treatments pl2]l3] we chose to 
Taylor expand this expression; here, we retain the full form to provide better accuracy in 
the case that 5'g*/g is not so small. 

The power spectrum (ll.2p is plotted in Fig. [3] for both positive and negative g. At 
sufficiently large momenta the spectrum rapidly becomes nearly scale invariant, with positive 
values of g resulting in a slight blue tilt and negative values of g yielding a slight red tilt. This 



/(^?eff) = /0(1 + ln(l/(/3^?eff)) + 0[gt^]), 



(2.8) 



(2.9) 
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behavior reflects the fact that the dual QFT becomes asymptotically free at high momenta, 
with the free theory itself corresponding to an exact Harrison-Zel'dovich spectrum. 

At lower momenta, the existence of the non-perturbative infrared scale giR becomes 
apparent, resulting in the peak/trough feature in the spectrum at g = egq^,. Note, however, 
that the perturbative calculation of f{g1^) underpinning the power spectrum ( 11.2^ breaks 
down when g^^ ~ gq^./q becomes of order unity (recalling that /i is a constant of order 
unity). This means that the perturbative result (11.21) becomes unreliable at low momenta 
close to the peak/trough feature in Fig. [3l Moreover, our approximation /3 = 1 is no longer 
justified in this regime and one should retain /3 as an independent parameter, see Section 
14. 3[ Since the smallest momentum scale appearing in the CMB is of the order 10"^ Mpc"\ if 
the power spectrum ( 11. 2p is to reliably fit the entire range of CMB scales, then we conclude 
that the maximum value of g is restricted to be of the order |(7|max ~ 2 x 10~^. We will use 
this rough estimate later in setting the prior for g. 

Note that the holographic model describes the very early universe; the end of this period 
should be the beginning of hot big bang cosmology. The asymptotic (late-time) metric 
obtained at the end of the holographic period satisfies Einstein's equations. This follows 
from the fact that the cosmologies we analyze are related to holographic RG flows via the 
cosmology /do main- wall correspondence [18], and the corresponding domain- wall spacetimes 
are known to satisfy Einstein's equations near their conformal boundary p^p5] . This then 
ensures the conservation of the standard gauge-invariant curvature perturbation which is 
the field dual to the trace of the stress energy tensor T. One thus expects a smooth transition 
to hot big bang cosmology, with the holographic period supplying the initial conditions for 
the subsequent evolution. 

Nevertheless, it would be very interesting to develop a detailed theory for the transition 
period, the analogue of the reheating period in conventional scenarios, and here we offer a 
few preliminary comments leaving a more detailed study for future work. In order to exit 
the holographic period we would need to modify the UV structure of the dual QFT (since 
the UV of the QFT corresponds to late times). This can be achieved by adding irrelevant 
operators to the QFT. At momenta far below the momentum scale quv set by the lowest 
dimension irrelevant operator, the computation of the 2-point function (and therefore of the 
power spectrum) is well approximated by the computation described above. Thus, as long 
as quv is much larger that the largest momentum scale seen by CMB, quv ^ 10~^ Mpc~\ 
the error made by omitting this period is very small. In principle, one could compute the 
corrections to the theoretical formula due to such irrelevant operators and extract from the 
data the best-fit value for quv- We leave such study for future work, but we note that the 
ability to fit the data well without these corrections indicates they are indeed smallH 

^Note however that the fit of the hofographic modef to data is better at low q than high see Section 
Perhaps the fit at high q would improve by including quv in the analysis. 
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3. Parameter estimation 



Let us now turn to a comparison of the model with the data. Holographic ACDM is described 
by the following free parameters, (f2Ji^, Qch"^, 6, r, A^, g). The first five are common with 
the usual parametrization of ACDMj; the sixth parameter is related to the coupling of the 
dual gauge theory, as discussed in the previous section. For now, we will take the power 
spectrum (11. 2p at face value. In the following section we compare the model to ACDM and 
calculate the Bayesian evidence after restricting g to values for which the assumptions used 
in deriving the power spectrum are satisfied. 

We perform our analysis for three data sets: the seven-year WMAP data [2T], and two 
combinations of data sets introduced in [21] as WMAP+BAO+i^o and WMAP+CMB. The 
former is a combination of WMAP7 with priors on the Hubble constant p2] and angular 
diameter distances [23]. The latter is a combination of WMAP7 with small-scale CMB 
experiments. As usual, to compute the posterior probability distribution on the parameter 
space of the model given the data, we use Bayes' theorem which, for a fiat prior probability 
distribution, relates it to the likelihood, i.e., the probability for the data given a choice of 
parameters of the model 

P{aM\D) ^ PiD\aM) = Ciau) , (3.1) 

where aM collectively denotes the parameters of the model M. For a given point in parameter 
space the likelihood can easily be evaluated, but the high dimensionality of the parameter 
space makes a grid-based approach prohibitive and other methods are required to sample it. 
We will present results from an analysis using CosmoMC [7f8] with the MultiNest sampler [9l 
[To] . We have also performed a more traditional Markov chain Monte Carlo analysis for the 
WMAP7 data set and found excellent agreement between the two. Both analyses made use 
of CAMB [2^f25] to evaluate the angular power spectra, version 1.5 of recfast [26] to compute 
the reionization history, and version 4.1 of the WMAP likelihood code p7j to evaluate the 
likelihood. 

As a consistency check, we first performed a parameter estimation for ACDM and found 
good agreement with the values determined by the WMAP team|l| Since we used a pivot 
scale = 0.05 Mpc~^, we repeat the parameters in Tabled] to facilitate comparison. 

The results for the mean values of the marginalized posterior distributions of the six 
parameters of the model and their uncertainties at 68% confidence level are summarized in 
Table [21 All of them are for a pivot scale = 0.05 Mpc~^, but notice that changing the 
pivot scale has only the trivial effect of rescaling the parameter g. Notice that for the best-fit 
values for g, the condition for parametric control gq*/q ^ 1 no longer holds for the lowest 

^Note the parameter is the ratio between the sound horizon at the time of last scattering and the 
angular diameter distance of the surface of last scattering. Physically, this ratio fixes the position of the 
acoustic peaks and 9 is tightly constrained by the data. Theoretically, 9 is a. function of ilbh^, ^Ich^ and h, 
so given fli,h^ and Q,ch? , it may be expressed in terms of the dimensionless Hubble parameter h. For more 
details see, e.g., Section 7.2 of [20] . 

^The slight but consistent increase in the value of the scalar spectral index can be traced to changes 
between recfast version 1.5 used in this work and recfast version 1.4 used in the WMAP7 analysis. 
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WMAP7 


WMAP±BAO±iJo 


WMAP±CMB 




0.02252 ± 0.00056 


0.02257 ±0.00053 


0.02265 ±0.00051 




0.1116 ±0.0054 


0.1127 ±0.0035 


0.1124 ±0.0048 


100^ 


1.0394 ±0.0027 


1.0400 ± 0.0026 


1.0411 ±0.0022 


r 


0.088 ±0.014 


0.088 ±0.014 


0.088 ±0.014 




(2.183 ±0.073) X 10-9 


(2.191 ±0.075) X 10-9 


(2.190 ±0.068) X 10-9 


Us 


0.969 ±0.014 


0.970 ±0.012 


0.969 ±0.013 



Table 1: This table shows a summary of the parameters of the ACDM model and their 
uncertainties at the 68% confidence level. 





WMAP7 


WMAP±BAO±/7o 


WMAP±CMB 




0.02310 ±0.00045 


0.02312 ± 0.00043 


0.02326 ± 0.00045 




0.1077 ±0.0051 


0.1120 ±0.0036 


0.1076 ±0.0042 


me 


1.0407 ±0.0026 


1.0406 ±0.0026 


1.0423 ±0.0022 


T 


0.087 ±0.015 


0.084 ±0.015 


0.088 ±0.016 




(2.146 ±0.088) X 10-9 


(2.172 ±0.086) X 10-9 


(2.151 ±0.084) X 10-9 


9 


-0.00127 ±0.00093 


-0.00136 ± 0.00094 


-0.00114 ±0.00088 



Table 2: This table shows a summary of the parameters of the holographic model and their 
uncertainties at the 68% confidence level. 

CMB momenta, and the data appear to push us to intermediate values of the coupling. We 
will return to this issue in Section 14.31 

The estimated values of the parameters shared between the holographic and conventional 
ACDM models all overlap in both scenarios. However, we see that the central value of Q^h"^ 
has moved by around one standard deviation. There are some degeneracies between the 
new parameters in the power spectrum and the parameters characterizing the background 
geometry. In Fig. |U we display these degeneracies in the form of a triangle plot for the 
WMAP7 analysis, and give the results for the individual variables in Fig. [51 

4. Model comparison 

We now turn to a comparison between the holographic model and ACDM. Two approaches 
have commonly been used to select models in the cosmology literature. The first is based 
on information criteria; these typically reward a larger best-fit likelihood while penalizing 
the introduction of extra parameters. The second approach, based on the Bayesian evi- 
dence, rewards the model that provides the best fit when averaged over its parameter space. 
Information criteria are relatively easy to evaluate since we already know the value of the 
likelihood for the best-fit point from the parameter estimation. The Bayesian evidence is 
harder to evaluate because it involves an integral of the likelihood function over the parame- 
ter space of the model, but it allows constraints on the parameters to be taken into account 



9 



0.12 
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0.08 
0.06 
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Figure 4: A triangle plot showing the parameters of the holographic model derived from 
WMAP7. The contours represent 68% and 95% confidence levels. 
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Figure 5: The marginalized likelihoods for the parameters of the holographic model as derived 
from WMAP7. 



more elegantly. This is especially important in the case of the holographic model since the 
form of the power spectrum (11.21) is only valid when the gauge theory is in the perturbative 
regime, hence there are natural priors that should be imposed on its parameters. We will 
discuss this in more detail in Section H^l we begin with a naive comparison using the best-fit 
likelihood. 



4.1. A first look: information criteria 

The most commonly used information criteria are the Akaike information criterion [28] 

A1C = -21n£bcst + 2A;, (4.1) 

where k is the number of parameters characterizing the model, and the Bayesian or Schwarz 
information criterion [29] 

BIC = -21n£best + A;lnA^, (4.2) 

where N is the number of data points. Since the holographic model in the regime of small 
coupling and the ACDM model possess the same number of parameters, both criteria reduce 
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Holographic Model 


ACDM 


A In £bcst 


WMAP7 


3735.5 


3734.3 


1.2 


WMAP+BAO+/Jo 


3737.3 


3735.7 


1.6 


WMAP+CMB 


3815.0 


3812.5 


2.5 



Table 3: This table summarizes the best-fit values for — In £ for both the holographic model 
and ACDM, as well as the difference between them. Positive numbers in the last column 
favor ACDM. The errors on our best-fit log likelihoods are estimated to be around 0.1. 



to a comparison of the value of the logarithm of the likelihood function at the best-fit point. 
We summarize the best-fit likelihoods for the two models in Table [31 

From Table IHl we may infer for example that, with all parameters tuned to their best-fit 
values, obtaining the observed WMAP data is approximately three times more probable 
according to ACDM than according to the holographic model. Note, however, that the 
probability of obtaining the data given the model with specific parameter values is not the 
same as the probability of the model given the data; it is the latter quantity that we really 
wish to evaluate and to which we now turn our attention. 



4.2. A closer look: Bayesian evidence 

Before presenting our results, let us quickly review Bayesian evidence. For a more detailed 
explanation, we refer the reader to [TT|12f30] and references therein. 

Our goal will be to compare the holographic model and ACDM, but we will keep our 
discussion general for now. A given model is labelled by M, and characterized by a set of 
parameters aM- In general, some of the parameters may overlap for different models, while 
others may be different. In the previous section, we determined the values aM for model M 
that are most likely given the data, i.e., maximizes P{aM\D). We now ask which model is 
most likely given the data, i.e., which choice of M maximizes P{M\D). Again we use Bayes' 
theorem to relate this to the probability for the data given the model, 

, , , P(D\M)P(M) , , 

P{M\D)= ^ ^^^^^ ' . (4.3) 

The unconditional probability for the data appearing in the denominator is a model-independent 
constant and may be eliminated in the comparison of two models by taking the ratio 

P{M,\D) _ P{D\M{) P{M{) 

P{M2\D) P{D\M2)P{M2) ■ ^ ' ' 

This is sometimes referred to as the posterior odds. The first factor on the right-hand side 
is often called the Bayes factor, while the second factor is the ratio of prior probabilities for 
the models under comparison. Typically, one assumes that all models are equally likely and 
this factor is set to unity. The quantity of interest is P{D\M), which is called the evidence 
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and often denoted E. The probability for the data given a certain model is of course nothing 
but the probability for the data given a certain choice of parameters, P{D\aM) = (^{oim), 
times the probability that this choice of parameters is realized, P(aA/), integrated over the 
space of parameters of the model. Consequently, 



Numerically, the evidence is difficult to evaluate because it is effectively a multidimensional 
integral over a function which is computationally nontrivial. However, an efficient solution 
to this problem was found by Skilling [31], and is implemented within MultiNest [9l[T0] . 

As in parameter estimation, we need a prior probability distribution for the parameters 
um- We will work with flat priors: namely, a prior probability that is constant over some 
defined region, and zero outside it. In this case, the evidence reduces to the integral 



where Vol a/ is the volume of the region in parameter space over which the prior probability 
distribution is non-zero. If we have a strongly peaked likelihood function which has support 
over only a relatively small region inside Vol^/, changing the prior region can strongly affect 
the computed evidence. Provided the changes to the overall volume of the parameter space 
do not add or exclude regions where C is large, the integral will be unaffected while Vol^f 
can change substantially and the computed evidence is inversely proportion to VoIm- 

With the exception of Hs and g both models have the same parameters. By using the 
same priors for the variables shared by the holographic and standard ACDM scenarios the 
ambiguity in the evidence associated with VoIa^ is minimized. However, the situation with 
g and Hs is more nuanced. For the holographic model we should restrict g to values where 
perturbative expansion used to derive (11. 2p is valid. There are several unknown factors of 
order unity in this perturbation expansion, so the permissible range of g is not specified 
precisely, but we therefore simply estimate the largest allowed value of the coupling Isflmax, 
permitting values of either sign. As discussed in Section [21 we determine the value of l^^l for 
which the naive expansion parameter \g\qif/q becomes unity at the longest scales to which the 
CMB is sensitive, yielding If^lmax ~ 2 x 10^'^. We then check whether the evidence depends 
strongly on this choice by setting \g\q*/q to values between 1 and 0.005 at the longest scales, 
the latter choice giving l^flmax ~ 10~^ and find that the evidence is only mildly dependent 
on the choice of If^lmax- 

The prior for Ug is less obvious since, unlike Is'lmax, it is a purely empirical parameter 
and we cannot restrict it by appealing to the internal consistency of some underlying theory. 
Moreover, our best information about Ug is derived from the WMAP data we are using to 
compute the evidence, and it would be inappropriately circular to set the prior on Ug directly 
from a parameter estimate derived from the WMAP data itself. On the other hand, choosing 
an overly generous range for Ug reduces the evidence and may overestimate the ability of 
competing models to explain the data. To account for this dilemma, we work with two 




(4.5) 




(4.6) 
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X 


^min 


^max 


Qbh'^ 


0.020 


0.025 




0.09 


1.25 


100^ 


1.03 


1.04 


r 


0.02 


0.15 




1.82 X 10^9 


2.71 X 10"^ 


Us (1) 


0.92 


1.0 


ris (2) 


0.9 


1.1 



Table 4: This table shows the priors used in the calculation of the evidence for the ACDM 
model. The two priors correspond to the narrow and broad choices discussed in the text. 

priors for n^, 0.92 < < 1.0 and 0.9 < Ug < 1.1. As is well known, the first choice includes 
essentially the entire range over which the likelihood is appreciably different from zero while 
the second is centered on the scale invariant Harrison-Zel'dovich spectrum. For these choices 
the ratio of VoIm is 0.08/0.2 = 0.4, shifting InE by about 0.9. In order to conclusively favor 
one model over another we would require a much larger difference in the respective values 
of the evidence, so this ambiguity will not qualitatively alter our results. For completeness, 
we summarize the priors on all the parameters of ACDM in Table 4. 

We also consider asymmetric priors in which the coupling g ranges from zero to a value of 
either +|5'|max or — |(7|max- We present results for \g\ma.x = 1 x 10""^ and also |(7|max = 2 x 10~^. 
The latter value is rather large, so we use it with the caveat that the assumptions implicit 
in our use of the holographic power spectrum may begin to break down for momenta that 
contribute to the lowest multipoles in the angular power spectrum. However, if we assume 
that the power spectrum can be trusted in this scenario we find moderate evidence suggesting 
g < using present data. Consequently, it seems worth exploring both the reliability of (11.21) 
at higher values of Is'lmax and attempting to determine whether fundamental theory predicts 
the sign of g, as answers to these questions could significantly improve our ability to constrain 
holographic ACDM, even with present-day data. 

We show the evidence computed using WMAP7 on its own and the combinations WMAP 
+BAO +Ho and WMAP+CMB in Figs. El [7] and El respectively. Along with the evidence 
for ACDM and the holographic model, the plots also show the evidence for the Harrison- 
Zel'dovich spectrum. As a guide to the eye, we shade the plots to indicate regions with 
Aln_E < 1 (white), 1 < Alni? < 2.5 (light gray), and 2.5 < A InE < 5 (darker gray). 
In the terminology of the Jeffreys scale, these intervals are "not worth more than a bare 
mention", "significant", and "strongly significant" [12]. On the basis of the narrower prior 
0.92 < < 1.0, we find that there is only weak to moderate evidence in favor of ACDM 
based on the data sets we have studied. 

When g is restricted to small values, the evidence we compute for holographic ACDM 
approaches the value found for the scale-free Harrison-Zel'dovich spectrum. This is easily 
understood: in this limit, the power spectrum (II. 2p is only weakly dependent on values of q 
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Figure 6: This figure shows the results of our calculation of the evidence for the WMAP7 
data set. ACDM is shown as a solid orange line and a dashed orange line for the narrow 
and broad priors on discussed in the text, respectively. The red line is the evidence for 
the pure Harrison-Zel'dovich spectrum. The green data points represent evidence computed 
for the holographic model, as a function of Ififlmax, as indicated on the horizontal axis. The 
red and blue points show the holographic model with asymmetric priors with positive and 
negative g, respectively. The shading indicates the difference in evidence according to the 
Jeffreys scale (see end of Section relative to ACDM with the narrow prior. 



which contribute to the CMB. Consequently, A^(g) depends very weakly on g, and dC/dg 
is close to zero. Thus the evidence integral effectively factorizes and the resulting J dg term 
in the numerator will cancel the corresponding term in VoIm- In this limit the expression 
for the evidence reduces precisely to the Harrison-Zel'dovich variant of ACDM. 

4.3. Sensitivity to the infrared scale 

As we discussed in Section |2l the 2-loop QFT calculation for f{g'^ff) is sensitive to an infrared 
scale giR ~ dyu^- This is reflected in the presence of the /2 term in (12. 6p . or equivalently 
the factor in (12. 8p . As long as one probes the theory on scales far above giR, for which g'^^ 
is very small, we have | lngq^/q\ ^ | In/sj. These terms therefore make only a very small 
contribution and their precise numerical values are of no significance. 

In our case, the scales that we probe are fixed: the WMAP momentum range is approx- 
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Figure 7: This figure shows the results of our calculation of the evidence for the 
WMAP-|-BAO+iJo data set, with the same conventions as Fig. |6l 



imately 10""^ ^ Q 10^^ Mpc^^ but the size of (^ym; equivalently g, is unknown. For 
sufficiently small g, the infrared scale giR is outside the WMAP momentum range and the 
theory is then insensitive to value of /s (since | \ngq^,/q\ ^ | In/sl for all modes in the CMB). 
This has been our reasoning in setting /s = |/i| thus far. 

For larger values of g, however, the infrared scale qi^ is inside the WMAP momentum 
range and the long wavelength modes in the CMB are sensitive to the infrared scale. In 
this case, it is not a good approximation to set /s = but we we should take /s = 
where the infrared parameter /3 encodes the location of the non-perturbative scale giR. In 
principle, /3 is fully determined by the dual QFT and is not an additional parameter of the 
theory. In practice, however, computing /3 is not an easy task and requires non-perturbative 
information {e.g., a large- iV resummation, see p^|TF] ). However, we can determine the value 
of the coupling for which the data becomes sensitive to (3. For very small values of g, the 
power spectrum (11. 2p used in the previous sections accurately reflects the predictions of the 
holographic model and is effectively independent of (3, so its precise values does not matter. 
However, there will be a range of values of the coupling for which both the spectrum is 
sensitive to the specific value of /3 and the perturbative expansion is still valid, and in this 
case our spectrum (II. 2p is incomplete. We determine the value at which /3 becomes important 
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Figure 8: This figure shows the results of our calculation of the evidence for the 
WMAP+CMB data set, with the same conventions as Fig. [6l 



by fitting to the power spectrum 

and computing the marginalized probability distribution for /3 with different values of If^lmax- 
The result is displayed in Fig. O As expected, for \g\max ~ 10~^, the marginalized 
probability distribution for /3 is approximately flat over several orders of magnitude, and /3 
is an irrelevant parameter. Raising If/lmax to 2 x 10~^, we see that the probability distribution 
for /3 is peaked near /3 ~ 1. For couplings \g\ of this magnitude, we are quite sensitive to 
the presence of the infrared scale. It is conceivable that the correct value is indeed of order 
unity, but it seems equally likely that /3 is a factor of ten larger or smaller than unity, in 
which case we are far from the result that would be returned by a precise calculation of the 
spectrum. The transition is not very sharp and occurs around \g\ ~ 10~^. We have seen that 
the best-flt points for the holographic model for the three data sets all lie at values of \g\ 
about an order of magnitude larger than this. It thus seems that the data may be pushing 
us into a regime of intermediate coupling where the predictions become infrared sensitive. 
This motivates further theoretical work to establish the predictions of the holographic model 
in this regime. 
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Figure 9: This plot shows the marginahzed probabihty distribution for the infrared parameter 
^ for various values of the coupling. The upper left corresponds to \g\max — 0.002, the upper 
right to \g\miix — 3.4 x 10""^, the lower left to jg'lmax = 5.8 x 10~^, and the lower right to 

l^lmax = 10"^ 
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5. Conclusions 



In this paper, we used the WMAP data to test the holographic ACDM model, and we 
conclude that the holographic model is compatible with current data. To compare the 
holographic model with conventional power-law ACDM, we evaluated Bayesian evidence for 
both models, paying close attention to the choice of priors. The prior on the coupling 
that controls the holographic model, is set to include only values for which the perturbation 
expansion underpinning the holographic power spectrum is valid over the entire range of 
momenta relevant to the CMB. For the power-law ACDM model, is a purely empirical 
parameter so the corresponding prior cannot be deduced on the basis of a purely physical 
argument. Consequently, we consider two different choices for prior; the first is a near optimal 
choice containing the entire range over which the likelihood is appreciably different than zero 
and the second is centered on the scale invariant Harrison-Zel'dovich spectrum. In both cases 
there is at most weak evidence for ACDM, relative to the holographic scenario, although the 
difference in evidence grows somewhat more pronounced if ground/balloon-based CMB data 
is added to the fit, or if we include measurements of baryon acoustic oscillations and Hq. 
More precisely, with the narrow prior of 0.92 < Ug < 1 for the ACDM spectral tilt the 
difference in — \n.E is of order 1.2 to 1.6, see Figs. El [7] and [HI while if one assumes the prior 
of 0.9 < Ug < 1.1, the difference in evidence is no longer significant. We conclude that at 
this point we do not find any strong evidence in favor of ACDM. 

We can also compare holographic ACDM to the exactly scale invariant Harrison-Zel'dovich 
power spectrum, and in this case the computed evidence for both models is essentially identi- 
cal, unless we allow the coupling g to approach values for which the underlying perturbation 
theory appears to become unreliable. Note that the holographic model with g = Q (i.e., 
free field theory in three dimensions) predicts an exactly scale invariant Harrison-Zel'dovich 
power spectrum, and is the only known microphysical model to do so, to our knowledge. 
Thus, the current data are also consistent with the dual theory being a free QFT. 

Above all, the results of this paper call for an improved theoretical understanding of 
the predictions of the holographic model. Our restriction on the maximum value of \g\, 
derived from a naive estimate of when the effective coupling g1^ becomes of order unity for 
the lowest CMB scales, appears to limit the holographic model to spectra with insufficient 
scale dependence. Besides accounting for the similar performance of the holographic model 
and the Harrison-Zel'dovich spectrum, this is consistent with the location of our best-fit 
estimate for \g\ near the top of the currently allowed range. A precise determination of the 
relationship between g and g1^, and hence of the maximum allowed value of \g\, is therefore 
important. Specifically, this relationship is given by fig^^ = gq*/q, motivating a full 2- 
loop QFT calculation of the constant f\. The result will in general depend on the Yukawa 
and quartic couplings, as well as the field content of the theory. Larger values for fi would 
permit a larger upper bound for the coupling \g\, allowing power spectra with a stronger scale 
dependence. For example, if /i ~ 10, then fig^g ~ 1 implies g'^^ ~ O(10~^), hence the three 
and higher-loop contributions are (relatively) small, even though the 2-loop contribution is 
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of order one. On the other hand, if g^^ is not small for all relevant momenta, the higher 
order terms would be needed exphcitly. 

Even if the 2- loop approximation is sufficient, there is another related issue that must 
be addressed. As we discussed in Section [2l the complete 2-loop results are sensitive to an 
infrared scale giR ~ gyu^^ leading to the /2 term in ( 12. 6p . or the /s factor in ( 12. Sp . On 
scales far above giR, these terms make only a very small contribution. The range of momenta 
that contribute to the CMB are fixed, so equivalently infrared effects are unimportant only 
when the coupling g is sufficiently small. To confirm this quantitatively, in Section 14.31 we 
parametrized the infrared effects via an undetermined constant /3. For \g\ less than ~ 10~^ 
the data are insensitive to the precise value of /3. However, this threshold is below the 
best-fit value l^^l ~ 10~^, so developing a better understanding of infrared effects is clearly 
important. In this regard, a direct computation of /3, either through large- iV resummation 
methods or through lattice simulations, appears the most promising line of enquiry. 

On the observational side, parameter estimates can be expected to tighten dramatically 
in the near future. In particular, the Planck satellite will measure the power spectrum 
accurately over a wide range of angular scales, with excellent signal to noise. Assuming 
broken scale-invariance at a similar level to the central values seen by WMAP, the Harrison- 
Zel'dovich limit of the holographic model will be ruled out with a very high degree of con- 
fidence. Empirically, WMAP measures the running in the usual running-ACDM model 
with roughly the same degree of precision as it measures the breaking of scale-invariance, 

— 1 [21]. The holographic model predicts that ~ jn^ — 1|, so it is perhaps unsurprising 
that the WMAP dataset does not permit us to make a strong distinction between the holo- 
graphic model and ACDM without running. However, Planck will measure a.^ at a much 
higher level of precision than the current central value of — 1, and future experiments 
will do even better [32j. Consequently, it appears that Planck will be able to rule out the 
holographic scenario discussed here if a running a.^ ~ — 1| is not detected [3]. In the event 
that such a running is detected, in order to distinguish the holographic model from running 
ACDM and other models permitting a strong running (for example, slow roll inflation with 
modulations [33]), it will then be necessary to examine higher logarithmic derivatives of the 
power spectrum. The ability of Planck to discriminate between such models is less clear and 
merits further investigation. 

One issue we have not addressed in detail here is the tensor power spectrum. The ratio 
r of tensor to scalar amplitudes predicted by the holographic model may be found in [Tp]: 
to leading order in g'^Q, r is a constant depending on the field content of the dual QFT. In 
particular, we emphasize that r is not parametrically suppressed, in contrast to models of 
slow roll infiation. We do not present detailed parameter fits including tensors, but using 
the WMAP dataset on its own with a tensor contribution proportional to A^(g) we find a 
weak preference for a nontrivial spectrum of primordial gravitational waves, consistent with 
the results of [13]. This is understandable, since tensor modes contribute to the tempera- 
ture anisotropies at large angular scales and combining these with a nearly scale invariant 
spectrum of scalar perturbations results in a better fit to the WMAP data. The estimated 
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distribution for r peaks at zero for a fit to the WMAP7+Ho+BAO dataset, so we do not 
attach any particular significance to this result. A primordial gravitational wave background 
is degenerate with and including this contribution drives the maximum likelihood value of 
g closer to zero. Moreover, Planck (and a variety of forthcoming ground and balloon-based 
CMB experiments) will break this degeneracy by tightly constraining the B-mode of the 
CMB polarization, in addition to the primordial temperature anisotropics. 

If future data provides strong evidence for holographic ACDM this would constitute the 
first observational evidence for holography, given that these models do not have a conven- 
tional realization in terms of weakly coupled gravity coupled to other fields. On the other 
hand, one should emphasize that the scenario discussed is a specific realization of holographic 
infiation: even the conventional ACDM model, realized in the usual way as an infiationary 
model, also has a holographic realization in terms of a strongly coupled QFT. Moreover, 
there are other holographic models based on deformations of conformal field theories (the 
first class in the discussion of Section 2), and these would in general have different obser- 
vational signatures. It would be interesting to extract the predictions for these models and 
confront them with observational data. In summary, these studies thus present a unique 
arena where theoretical ideas about Planck scale physics can be tested observationally: an 
exciting period lies ahead of us. 
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